Quantum squeezing and entanglement in a two-mode Bose-Einstein condensate with 

time-dependent Josephson-like coupling 
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Dynamical evolution of quantum mechanical squeezing and entanglement in a two-mode Bose- 
Einstein condensate (TBEC) with an adiabatic, time- varying Raman coupling is studied by finding 
analytical expressions for these quantities. In particular, we study the entanglement between the 
atoms in the condensate as well as that between the two modes. The nature of the enhanced 
quantum correlations in TBEC is clarified by considering squeezing and entanglement both with 
and without the nonlinear interaction turned on; it is found that entanglement approaching maximal 
value can be achieved with the nonlinear interactions present. Somewhat counter-intuitively, greater 
squeezing is found in the absence of nonlinear interactions. This is due to the collapses and revivals 
of the TBEC quantum state induced by the nonlinear interactions. In addition, results involving the 
self-trapping phase state of TBEC indicate potential for creating a dynamically stable, macroscopic 
entangled quantum state which is relatively robust with respect to atom number fluctuations. 

PACS numbers: 03.75.Kk,03.75.Mn,03.75.Gg 



I. INTRODUCTION 

Bose-Einstein condensates (BECs) provide a useful 
system for investigating matter wave quantum squeez- 
ing and entanglement. Squeezing in BECs has already 
been investigated previously from a number of different 
perspectives [1|, 

H HH Squeezed states are quantum 
states for which no classical analog exists 0. The defini- 
tion of quantum mechanical squeezing, or the reduction 
of quantum fluctuations below the standard quantum 
limit (SQL), derives directly from the Heisenberg Uncer- 
tainty Principle: For two arbitrary operators A and B 
which obey the commutation relation [A, B] = C, quan- 
tum squeezing exists when one of the variables satisfies 
the relatio 
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[Ai] 2 <-|(C)| 



(1) 



Although the origin of quantum squeezing in many ex- 
periments is enhanced quantum correlations, it must be 
noted that the inequality JQ) is not a sufficient proof of 
enhanced quantum correlations between the particles, as 
discussed by Kitagawa and Ueda in the context of spin 
squeezing |8|. A measure of the amount of spin squeezing 
was defined as|l,|3 



2[AJi 



(2) 



where [A J±] 2 is the variance in the direction orthogonal 
to the total spin vector for an iV spin-i particle system. 
This definition takes into account the effect of quantum 
correlations as it derives from the fact that a minimum 
uncertainty state for TV elementary spin-i particles has 
spin J/2 (J = N/2) equally distributed over any two 
orthogonal components normal to the mean spin direc- 
tion. This follows from the fact that, in the absence of 
quantum correlations, the total variance in the normal 
direction is simply given by the sum of the variances of 



the individual elementary spins. Any state with uncer- 
tainty less than the minimum uncertainty state, i.e. spin 
squeezed due to enhanced quantum correlations between 
the particles, therefore gives £ < 1. 

Quantum entanglement, which is closely related to 
quantum squeezing, has recently generated a lot of ac- 
tivity amongst researchers, owing to its significant role 
in the studies of quantum information theory. A large 
number of theoretical studies exist on generating entan- 
gled states using BECs in a variety of different phys- 
ical settings. Many of these involve two-mode BECs 
(TBECs) in which atoms in two different hyperfine states 
are entangled [Tfl ITU Il2l IT3L IT^ | . In simple terms, a quan- 
tum state is said to be entangled when the state cannot 
be written as a simple product state, the most promi- 
nent example being the famous Bell-state for a two par- 
ticle system: \ip) = (|1)|0) + |0)|l))/-y/2, where "1" and 
"0" represent spin "up" and "down" states respectively. 
From the inseparability criterion for the A-particle den- 
sity matrix, a parameter for the enta ngle ment between 
atoms in TBECs has been derived [HJ E3, 
form similar to Eq. J2J): 



which has a 



A [Am • Jf 



(n 2 • J)2 + (n 3 • J) 2 



(3) 



where m, i = 1, 2, 3 are unit orthogonal vectors, and J is 
the well-known Schwinger angular momentum operator 
representing a many particle system. It is easy to see 
that Eq. © 

may be viewed as a generalization of Eq. 

On the other hand, it has also been argued that quan- 
tum entanglement in a TBEC is more meaningful when 
one considers the system as a bipartite system of two 
modes, analogous to considering entanglement of electro- 
magnetic modes, as opposed to the photons themselves, 
as there is no definitive measure for entanglement be- 
tween three or more subsystems |l 3 ■ In addition, the 
two modes of a TBEC are clearly distinguishable and 
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experimentally accessible subsystems. In this case, the 
standard measure of entanglement is the von Neumann 
entropy of the reduced density operator of either of the 
subsystems: 



E(pa) = -Trp A log 2 (pA), 



(4) 



where pA is the reduced density operator for subsys- 
tem A defined by the partial trace over the other sub- 
system B, pa = Trap. It can be easily shown that 
for a general bipartite system of TBECs written in the 
Fock basis \ip(t)) = J2 m c m{t)W ~ ™>) A \m) B , Pa = 



JY 



|to)(to| and hence the dynamically evolv- 



ing bipartite entanglement in TBECs can be parameter- 
ized by 



E(t) =-J2 |c m (t)| 2 log 2 |c m (i)| ; 



(5) 



In this paper, we calculate analytically dynamical evo- 
lution of quantum squeezing as well as quantum entangle- 
ment for both types of system decompositions of a TBEC: 
N particle subsystems based on Eq. Q and the subsys- 
tems composed of the two modes based on Eq. (JSJ). We 
shall consider a TBEC with a continuous time-dependent 
Raman coupling between the two states, with and with- 
out nonlinear interactions. A physical realization of such 
a system already exists; the Josephson-like coupling be- 
tween the two levels is provided by an external laser in 
the Rb TBEC of Ref. [16j , while a state-dependent mag- 
netic field gradient may be applied to induce Josephson 
tunneling in the Na spinor system of Ref. • We con- 
sider a system with an adiabatically time-varying cou- 
pling with an off-resonance atom-light interaction. The 
TBEC is itself an interesting and important system to 
study: not only does it display effects similar to that of 
Josephson junction in superconductors, such as the col- 
lapse and revival phenomenon in macroscopic scale, it is 
also the simplest Bose-Hubbard model with a two site 
lattice potential for which complete description of entan- 
glement can be given. 

The central result of this paper is that we use an ex- 
act solution to the time-dependent Schrodinger equation 
to find analytical expressions for the dynamical evolu- 
tion of squeezing and entanglement. In particular, we 
explicitly identify the effect of the nonlinear interactions 
on squeezing and entanglement. Since one of the com- 
ments of Ref. 15] with regard to the atom-atom en- 
tanglement parameterized by Eq. © was that of prac- 
ticality, we shall restrict our calculations to two differ- 
ent possible initial states for a TBEC, the Dicke state 
and the phase state. We consider atom-atom entangle- 
ment for which clear physical measurement is in principle 
possible, namely those involving the measurement of the 
mean and variance of the atom number difference and 
the relative phase between the two species. 

The paper is organized as follows: In Sec. II, we de- 
scribe, using the Schwinger notation for the angular mo- 
mentum operators, the TBEC and the exact solution to 



the time-dependent Schrodinger Equation. In Sec. Ill, 
we present our main results: the dynamics of the TBEC 
system with a time-dependent Josephson-like coupling 
between the two species and the evolution of quantum 
squeezing and entanglement for the two types of subsys- 
tem decompositions for various values of laser couplings, 
with and without the nonlinear interaction. We conclude 
in Sec. IV. 



II. FORMALISM 

A. Hamiltonian 

We consider atomic BECs in two different hyperfinc 
states trapped in a single trap, with a time-varying Ra- 
man coupling between the two levels given by a spatially 
uniform electromagnetic field. The Hamiltonian for this 
system, with the annihilation operators for the two dis- 
tinct states denoted a and b under the two mode approx- 
imation is the following: 



H — H a + Hb + H int + -Hlas 
H a 

H, 

Hint 

Hiss 



. L/ a a + - 1 * - 

uj a a ] a H — —a' a 1 aa 
^a+aStS 



(6) 
(7) 

(8) 

(9) 
(10) 



H a and Hb describe the two condensates undergoing 
self-interaction while Hint and Hi as describe the conden- 
sates interacting with each other via collisions and laser- 
induced interactions respectively. _ffi as describes a time- 
dependent coupling, rendering an overall time-dependent 
Hamiltonian for the system. The Hamiltonian may be 
rewritten by employing the Schwinger notation for the 
angular momentum operators, namely, J x — ^(a^b+b'a), 
Jy = hic^b — b*a), and J z = ^(a)a — wb) with the 

Casimir invariant J 2 = =j(t + 1) where N = a) a + wb 
is the total number operator and is a conserved quantity. 
Physically, eigenvalues of the operator J z represents the 
difference in the number of atoms in different hyperfine 
levels, while J x and J y takes on the meaning of relative 
phase between the two species. The z component an- 
gular momentum eigenstates, known also as the Dicke 
states [Tsl llflj . are written as \j, m) where m = —j, . . . , j 
with J z \j, m) = m\j,m). Here, j = N/2 and is the quan- 
tum number of angular momentum. The raising and low- 
ering operators are defined in the usual way as J± = J x ± 
iJ y such that J±\j, m) = \J (j =F to) (j ± to + 1) \j, mil). 
In terms of the angular momentum operators, the Hamil- 
tonian takes the form 20, 21]: 



H (t) = ljq J z + qJl + Sl(t) J+e 



(11) 
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where uj = u a - uj h + (N- \){U a - U b )/2, q = (U a + U b - 
U a b)/2. It is noted that U a (Ub) or U a t may, in principle, 
be tuned via Feshbach resonance through the application 
of an external magnetic field|22|: the factors u>o and q are 
consequently adjustable parameters. 



B. Solution to the time-dependent Schrodinger 
Equation 

An exact solution to the Schrodinger Equation 



ih-m)) = m)\m) 



(12) 



with the time-dependent Hamiltonian H (t) of Eq. 

can be given in terms of a time evolution operator U (t) , 



U{t) = ft(t)e- iH 'WR(0), 



(13) 



where R is a time dependent unitary transformation de- 
fined as 



R(t) = exp 



and 



H'{t) = RH{t)R} - iR—R). 



(14) 



(15) 



The parameter A in Eq. I|14|) is an auxiliary parame- 
ter which may be chosen to simplify the transformed 
Hamiltonian, H' . It can be shown that R generates a 
gauge transformation under which the time-dependent 
Schrodinger Equation is covariant 23J : 



of the system. We note here that by ignoring the 
two-photon transitions for the purpose of obtaining 
analytical solutions, the effect of "two-axis counter- 
twisting" Q is not included in our effective Hamiltonian. 
For concreteness, we shall consider in this paper the case 
ip(t) = At where A gives the detuning of the laser from 
the \A) — > \B) transition between the two species. Also, 
we shall consider the case U a — Ub > U a t and identical 
trapping potentials for the two species, corresponding to 
the \F = 1, Mp = ±1) hyperfine states of Na trapped in 
an optical dipole trap [l(| • This implies q > as well as 
lu = i.e. A = -7r/2 in Eq. [fTij) . 



C. Initial quantum state of TBEC 

We consider in this paper an SU(2) atomic coher- 
ent state or a coherent spin state (CSS), \9, </>)0, 
as the initial quantum state for a TBEC. A CSS is a 
minimum uncertainty state that describes a system with 
well-defined relative phase between the two species, and 
provides a good description of TBECs under suitable ex- 
perimental conditions [25j. It is defined mathematically 
by applying the rotation operator on the extreme Dicke 
state \j,j) or The definition that we use in this 

paper is: 



exp 



2 (J - £ 



(19) 



A CSS \0,4>) is therefore an eigenstate of the spin com- 
ponent in the (0,<j)) direction Jg^ = J x sin 9 cos <j> + 
J y sin 6 sin 4> + J z cos <fi with eigenvalue J where 9 and 
4> denote polar and azimuthal angles. It can be shown, 
following from Eq. H19J) . that CSS may be written as a 
superposition state: 



ih-W{t))=H'{t)W{t)), 



(16) 



where the transformed state \ip'(t)} = R\ip(t)). For our 
case, A is chosen such that the Hamiltonian is diagonal 
in the J z representation. With an additional assumption 
of adiabaticity conditions dip/dt ~ and dX/dt « 0, 
and also assuming that the two-photon transition terms 
proportional to J+e l2(p ^ and J? e - l2i p{t) can b e neglected 
l24l. one obtains for the transformed Hamiltonian 



H'(t) 



where A is chosen such that 
tan A = - 



(17) 



(18) 



The time evolution of any observable A is then 
given in the Heisenberg picture as (A(t)) = 
(ijj(0)\ W (t)Au (t)\iP(0)) , where U(t) is defined in 
Eqs. (|13I15|) and 1^(0)) is the initial quantum state 



E 



TV 



(20) 



where 



,{9 A) = (Cf| m )^cos^Q)sin^^ 



xe 



i(j-m)cj> 



(21) 



and C™ denotes the combination, C,™ = n\/\{n — m)!rn!]. 
Further extensive discussions on the properties of the 
CSS may be found in Refs. 0,0]. In this paper, we 
shall consider initial CSS's that may be produced exper- 
imentally, \9 = 0,0 = 0), and \9 = tt/2,0 = 0). The 
CSS |0 = 0, (f> = 0) is simply a Dicke state \ j,j) in which 
all atoms are found in one hyperfine level, and the CSS 
1 9 = 7r/2, cf> = 0) is an eigenstate of J x which is a state 
with a well-defined phase difference 4> — 0, and is a phase 
state of a two-mode boson system 26\. States with differ- 
ent values of 9, although simple to deal with mathemati- 
cally, are difficult to produce and observe experimentally. 
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For the same reason, we shall restrict our consideration 
of squeezing to those in the x, y, and z directions since 
the variances in the relative phase and atom number are 
physical quantities clearly measurable without the am- 
biguity associated with the measurement of variances in 
other directions. 



III. RESULTS 



based on Eq. JT3J, 

(j a ) = (d,<t)\R\0)e iA 'R(t)J a &{t)e- i6 'R(0)\6,<p). 

In the first instance, this may be calculated by using the 
commutation relation amongst the spin-J components, 
namely [J z , J±] — ± J± and [J + , J ] = 2J Z , and employ- 
ing the Baker-Hausdorff theorem[l9|. Useful identities 
are given in Eqs. I|23I25|I with ip(t) = At: 



In this section, we present various results for the two 
initial states for a TBEC, the Dicke state and the phase 
state. The coefficient q provides the strength of the scat- 
tering interaction between bosons, and its magnitude in 
relation to the tunnelling coupling (as we consider 
luq = in this paper) is an important parameter in de- 
termining distinct coupling regimes |27j. Of the three 
coupling regimes, namely, the Rabi (q/fl <C 1/-ZV), the 
Josephson (1/N < q/fl < N) and the Fock (q/Q > N) 
regimes, it is clear that the regimes of more physical in- 
terest are the Rabi and the Josephson regimes, since, in 
the Fock regime the tunnelling coupling is overwhelmed 
by the nonlinear interaction term, resulting in no new 
effects of particular interest. This has been confirmed 
numerically |23| • 

In the Rabi regime, the tunnelling coupling dominates 
with g/fi near zero. We shall therefore consider the lim- 
iting case of q = i.e. the linear regime which not only 
illustrates the essential physics of the Rabi regime but 
also demonstrate clearly the effect of nonlinear interac- 
tion on squeezing and entanglement. In addition, we shall 
consider the lower end of the Josephson regime, as it 
was found that important features rapidly disappear as 
one nears the Fock regime. The two couplings consid- 
ered in this paper are therefore q = and q = qj with 
qj/fl = 3/N for the Josephson regime, where we choose 
N to be 400. We also have an additional parameter in 
our model not considered in some of the previous work, 
that of a detuning A. As will be shown below, this en- 
ters mainly as a phase shift which affects the results in a 
nontrivial manner. To illustrate the effect of this param- 
eter we shall consider detuning of the same magnitude 
as the q for the Josephson regime, (A = qj) as well as 
a relatively large value of detuning (A/J7 = 10g 3 ). In 
particular, the results in the linear (q — 0) regime will 
help us see the role of A with respect to q. The result for 
A = can be easily interpolated from these two values 
of detunings so will not be explicitly presented. 



A. Dynamics of the macroscopic spin vector 

We first consider how the quantum state of a TBEC 
represented by the fictitious spin vector evolves in time 
by calculating the expectation values (J a (t)) where a 
stands for x, y, and z. This later gives us insight into how 
the squeezing and entanglement in a TBEC dynamically 
evolves. We evaluate the mean values of spin components 



■At 



-2iAt\ 



RJ z R f = -^(J+e 

Rj+tf = hj+-L 

RJ-R) = i(J_ - J+e 2lAt ) + j : 



J_e- iAt ), (23) 



J z e 

Al 



-iAt 



(24) 
(25) 



In principle, it is possible to repeatedly apply the Baker- 
Hausdorff theorem to obtain a general expression for 
J a (t) = U^(t)J a U(t); for the special case of q = 0, the 
expressions may, in fact, be simplified further, giving 



and 



N 

(J z (t)) = — cos 9 cos[20i + At 



N 



(26) 



(4(*)> = Y e T!At 



^fi sin 



cos ( 


9-tt/2\ 




2 J 



' — tt/2 



sin(2fii + At 



(.27) 



It is found however that the expressions quickly become 
unwieldy especially for q ^ 0. The way around the prob- 
lem is to note an important and useful point that, for the 
Hamiltonian under consideration, the transformation op- 
erator R happens to have the same form as the generator 
of SU(2) atomic coherent states[Eq. I|19|) ] if one writes 
ip(t) = —(f>. In particular, the effect of R is to rotate the 
CSS such that R(0', <f/)\0, <f>) -»■ \6 + 6', <j> - 0'), i.e. 



R(t)\0,4>)= T*L{8-*l*,<l>-*t)\3,m)> (28) 



m=-j 



where 7t 3 m (9, 



is defined in Eq. I|21|) . Using the well- 
known relation for the raising and lowering operators J± , 
and the fact that F(J z )\j,m) = F(m)\j,m) where F de- 
notes some analytic function, one can obtain, after some 
algebra, an analytical expression for (J z (t)} as a summa- 
tion: 

N/2-1 

(lit)) = - V^m) 
m=-JV/2 



X COS 



2Qt - q ( m + - )t + At 



,(29) 
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A = q, 
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FIG. 1: (J z (t))/N, (J x (i))/N, and {J y {t))/N for q = 0. Solid 
line: initial Dicke state, \8 — 0, cj> = 0} Dashed line: initial 
phase state \8 — n/2,<f> = 0). Left column: A = qj. Right 
column: A = Wqj , where Qj/Q = 3/N 



FIG. 2: (J z (t))/N, (J x (t))/N, and (J y (t))/Nior q = qj. Iden- 
tical details as in Fig. 



where we have defined 



V a (9,m) = C 



N 

N/2+m+l 



t + m + 1 



X cos 



2iV 



-tt/2 



x tan 



N—2m—a 



e-Tr/2 



(30) 



Following somewhat more involved but identical steps as 
above, we obtain: 

N/2-1 

(J±(t)) = e Tl(A *-" /2) Pi(0,m) 

m=-N/2 



x sm 



29.1 - <i { m+ - ) l + A/ ■ 



N/2 

e* iAt V (6, 

=-N/2 



j - m 



(31) 



from which (J x (t)} and (J v (t)) may be deduced: J x — 
i(J+ + J_) and J a = ^(J + - J_). 

We first present in Figs. ^ and [21 the expectation val- 
ues (J a (t)), a = x,y,z scaled to the total number of 
atoms N for q = 0, and q = qj respectively for the de- 
tunings A = qj and A = lOqj. It is noted that no sig- 
nificant qualitative difference is observed in these scaled 
amplitudes when one uses different number of atoms. In 
all the figures in this paper, the solid line represents the 
result for the initial CSS \0, = 0) with 9 = i.e. the 
Dicke state, while the dashed line corresponds to the ini- 
tial phase state, 9 = tt/2. In Fig. ^ it is seen that, for 



the initial state \6 = 0, <f> = 0), the y and z components of 
spin, (J y ) and (J z ), undergo oscillatory evolution with a 
tt/2 phase shift. This implies that the macroscopic spin 
vector may be visualized as if undergoing a circular mo- 
tion in the y-z plane. Superposed on this motion is the 
gradual increase of the oscillation in the x component, 
( J x ) with the amplitude of the y component being pro- 
portionally reduced with a tt/2 phase shift. This implies 
that the spin vector initially undergoing a circular motion 
in the y-z plane rotates around the z-axis with increas- 
ing amplitude. This motion gets more pronounced for the 
larger A in that the spin vector undergoing circular mo- 
tion in the y-z plane rotates around the z-axis relatively 
quickly into a circular motion in the x-z plane. This 
continues on to return to the y-z plane, repeating this 
pattern over time. For the initial state \9 — tt/2, f/> = 0), 
there is no oscillation in the y-z plane, and only a slow 
rotation through the x-y plane is observed. This state 
may be identified as a "self-trapping" state, as there is no 
transfer of populations during the evolution. With higher 
A, the frequency of this rotation is clearly increased. 

For the case of q = qj, Fig. [3 shows the dephasing or 
"collapse" of oscillations for the initial state \0 = 0, <fi = 
0). The three spin components collapse to give average 
spin of almost zero. As will be seen below, this has a 
significant impact on the variance, and consequently on 
squeezing and entanglement in a TBEC. For the self- 
trapping initial state \9 = 7r/2, <f) = 0), there is no change 
in the behavior with q = qj. As to be expected, higher 
values of nonlinearity q results in shorter time scales for 
the collapse and with higher number of atoms, no sig- 
nificant qualitative difference is observed except that the 
collapse happens somewhat faster. The long time simu- 
lation is provided in Fig. |3 demonstrating clear revivals. 
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tive squeezing due to the rotation of the coordinate axes, 
while squeezing observed for q = qj comes from a com- 
plex combination of coordinate rotation on top of the 
interatomic collision effect of the one-axis twisting term 
proportional to J 2 |8| • In order to calculate the amount of 
squeezing in the system, we first calculate the variances 
in the three components of the macroscopic spin vector, 
i.e. [AJ Q (i)] 2 = - (J Q (t)) 2 , a = x,y,z which, 

as mentioned above are experimentally accessible quan- 
tities. The unitary property of R i.e. WR = 1 is used 
to evaluate terms of the form RJ a JpW = //•/,, If h'J : h" 
along with the identities Eqs. (|23I25|) . 
It can be shown that (J^(t)) is given by 



200 400 
nt (Jtrad) 



200 400 
£2t (icrad) 



FIG. 3: (J z (t))/N, (J x (t))/N, and (J v (t))/N for q = qj, sim- 
ulated for a longer time to observe revivals. Identical details 
as in Fig. Q 



The zero mean spin due to collapse indicates that the 
quantum state of a TBEC evolves from a CSS which is 
similar to the usual coherent state in quantum optics to 
a state which is similar to the number eigenstate with an 
equal number of atoms in each mode. The state that the 
CSS evolves into is, however, not exactly a self-trapping 
state as it clearly revives back into a CSS; it may be 
viewed as a quantum state with a small spread around the 
mean self-trapping state. For lack of better terminology 
we shall refer to this collapsed quantum state as a " quasi- 
self-trapping state" in this paper. 



B. Quantum squeezing and fluctuations 

In the language of spin squeezing, any reduction of 
the quantum noise for q = corresponds to an effec- 



1 N/2 ~ 2 f N 

(J 2 At)) = g E W,m)( T -m-l 

m=-iV/2 

x cos[4Q£ - 2q(m + l)t + 2At + 20] 

N/2-1 



+7 E V 2 (0, 



m=-N/2 



N 

m ) | m 

2 



N 



+V (e,m) ( — + m + l 



(32) 



so that the variance in the z direction, [AJ z (t)] 2 , is sim- 
ply given by subtracting the square of Eq. (|29|l from Eq. 

On the other hand, (J 2 ) = \[(j 2 + ) + (J+ ,L) + 

(J-J+)] and <J y 2 ) = |[<J + J_) + (J- J+) - (J 2 ) -(J*)]. 

The required expectation values of ( J±) and (J±J T ) are 
given by: 



N ^ 2 e -2,At , N 
{J+(t)) = 22 — o — v 2(0,m) I — - to - I ) cos[40t- 2q(m+ l)t + 2At + 20] 



-—N/2 
V/2-] 

E 



N/2 1 e -2iAt ( (N \ (N 

—— i V 2 (9, TO) --to + V Q (6, to) I — + TO + 1 
-N/2 IV/ V 



+ (2m + l)i sin 



21 >/ -<l[m+^)l + A/ - 



N/2 

+ Yl e~ 2lAt Vo(0,m) 

J m=-N/2 



J — TO 



(33) 
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(J+(t)J-(t)) 



N 



N/2-2 ^ 

^2 7p2(®, m) I y - m - 1 ) cos[40i - 2q(m + l)t + 2 At + 2<f\ 

-n=-N/2 ^ 



m=-N/2 v 

42?i (6»,m) cos 



2fii -g(m+-)i + At + 



iV/2 
m=-JV/2 



j - m 



(34) 



with (J 2 (t)) = (J+(t))* where asterisk (*) denotes 
complex conjugate, and (J_(t)J+(t)) = (J + (i)J_(i)) — 
2(J z (t)}, as to be expected from the commutation rela- 
tion [J_|_, J_] = 2J Z . From these expressions and those 
for (J x (t)) and (J y (t)) obtained above, the variances 
[AJ x (t)) 2 and [AJ y (t)] 2 may be calculated. 

We find that, for q = 0, the variances [AJ a ] 2 , a = 
x,y,z, show an oscillatory behavior which is bounded 
above by N/4 [Fig. 0]. This is an expected result for a 
CSS, which is known to have variances in the standard 
quantum limit (SQL) of J/2. One key observation re- 
garding Fig. 0]is that the variances do go below the SQL. 
This kind of squeezing which is due to the rotation of the 
coordinate axisjg naturally occurs since a spin vector is 
an eigenstate of spin in one direction with zero variance 
in that direction. As the spin vector traverses the phase 
space due to dynamical evolution, the error ellipsoid fol- 
lows the path of the spin vector in such a way that the 
minor axis of the ellipsoid periodically lines up with the 
x, y or z axis, resulting in the reduction of quantum fluc- 
tuations in that direction. For small A, the variance in 
the z direction show oscillatory reduction with 7r/2 phase 
shift with that of the y component. The variance in the 
x direction is seen as more or less maintaining its value 
near the SQL. For larger A, the uncertainty in the x 
direction is found oscillating with increasing amplitude, 
which is consistent with the error ellipsoid following the 
motion of the rotating spin vector discussed above. 

For the q = qj case, the variances in the x, y, and z 
directions give significantly different behavior, as shown 
in Fig. [3J The unexpected feature is the large variance 
equivalent to the maximum relative uncertainty (stan- 
dard deviation) of order ±35% per measurement. Math- 
ematically, this can be understood as the corollary of 
the collapsing mean spin components; the collapsing spin 
vector implies that the the variance [A J a ] 2 = (</ 2 ) — ( </ Q ) 2 
is dominated by the (J 2 ) term of the order TV 2 as it can- 
not be cancelled by the ( J a ) 2 term. Careful analysis of 
the variance for the q — case discussed above reveals 
that the terms of the form (J 2 ) are almost exactly can- 
celled by the ( J a } 2 which are also of the order TV 2 to 
maintain the variances of the order of the SQL, N/4. 
This type of cancellation cannot happen as the spin vec- 
tor collapses with q = qj. Physically, the increased un- 
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FIG. 4: Variances \AJ z (t)} 2 , [AJ x (t)] 2 , and [AJ y (t)] 2 for q = 
0, scaled to SQL of N/4 for clarity. Scaled variance of 1 in 
this figure corresponds to the SQL. Identical details as in Fig. 
□ 



certainty can be attributed to the quantum state evolv- 
ing away from the minimum uncertainty state of CSS 
into a quasi-self-trapping state with much higher quan- 
tum fluctuations. In order to display periodic behavior 
in the variances, we plot in Fig. EJthe variances over the 
identical time period as in Fig. [3] The result is in agree- 
ment with the purely numerical result of Tonel et al. |28( . 
It is found that, as the revivals occur, the variances do 
go below the SQL periodically. Although it is not clearly 
depicted due to scaling, the results for the initial phase 
state \9 = tt/2,(/) = 0) are identical to that for the q = 
case i.e. the variance stays at 1 for the atom number 
difference, and oscillates between and the SQL for the 
relative phase. 
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FIG. 5: Variances [AJ z {t)\ 2 , [AJ x (t)] 2 , and [AJ y (t)] 2 for q = 
qj, again scaled to the SQL, JV/4. Identical details as in Fig. 
□ 
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FIG. 6: Variances [AJ z (t)] 2 , [AJ x (t)] 2 , and [AJ y {t)] 2 for q = 
qj scaled to the SQL, N/4, simulated for a longer time to 
observe revivals. Identical details as in Fig. Q 



C. Dynamical evolution of quantum entanglement 
in a TBEC 



As mentioned in Ref. the amount of entangle- 

ment depends on the way a system is partitioned into 
subsystems. In this subsection we calculate the dynam- 
ical evolution of quantum entanglement in two different 
subsystem decomposition of the TBEC: entanglement be- 
tween the particles and that between the two modes. 



1. Entanglement between the particles 

The atom-atom entanglement is parameterized by Eq. 
©. Equation J2Jl is useful since, regardless of the actual 
spin direction at a given time, one need only to identify 
three unit orthogonal directions to check whether the 
system has become entangled due to system dynamics. 
If one's goal is to identify maximum entanglement possi- 
ble between the atoms, it is necessary to scan through all 
possible rii's within the unit sphere at each point in time. 
However, it should be noted that keeping to a fixed di- 
rection may help in constructing an actual experimental 
scheme to measure ^{t). In particular, we consider in 
this paper experimentally meaningful variances and am- 
plitudes of the macroscopic spin vector in the x, y, and 
z directions. We therefore consider 



ejt) 



N[AJ a (t)}' 



(Mt)) 2 + (MtW 



(35) 



where a, (3, 7 cycle through x,y, z. 

For q = 0, it was found that = 1, a = x, y, z, at 

all times for both initial states i.e. not entangled despite 
the fact that, as already seen in Fig. squeezing due 
to the rotation of the coordinate axis does occur. This 
shows that the observed squeezing is not accompanied 
by an increase in quantum correlations between the par- 
ticles. On the other hand, for q = qj, it is found that ^ 
is maintained at 1 for the initial phase state as to be ex- 
pected, while for the initial Dickc state, the very quickly 
increasing quantum fluctuations and the collapsing spin 
vector due to dephasing £^(t) large. Physically, this is 
due to the fact that the CSS is evolving into a quasi- 
self-trapping quantum state with low atom-atom corre- 
lations. A correct balance between the rapidly increasing 
variance and the collapsing spin vector is clearly needed 
in order to maintain entanglement with < 1. From 
Fig. [SI it is seen that the variance in J x increases most 
slowly amongst the variances, while Fig. [5] indicates the 
spin vector ( J y ) + (J z ) oscillating with a non-zero ampli- 
tude, showing the most promise in finding entanglement 
in this spin direction. In Fig. \7\we show £jl(t) for q = qj, 
for A = qj and A = 10<7j. Differently from other results, 
it was found that with A = 0, > 1 for all times, in- 
dicating the crucial role A plays in this system. For the 
initial Dicke state the parameter quickly becomes large, 
clearly indicating that quantum correlations between the 
atoms are destroyed rapidly. However, it is seen that the 
parameter £jl (t) does dip below 1 for a brief time period 
demonstrating atom-atom entanglement. 

Somewhat counter-intuitively, one may extend the du- 
ration over which the atoms are entangled by reducing the 
nonlinearity constant q; it is clear that with a smaller q, 
the variance in J x increases more gently while the col- 
lapse of the spin vector occurs over a longer time scale, 
extending the duration over which f£ < 1. This is pre- 
sented in Fig. |HJ for (7 ~ qj /10- As before, the increase 
in A results in the increase in oscillation frequency. This 
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FIG. 7: Entanglement parameter £ x for q — qj. Identical 
details as in Fig. 
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FIG. 8: Entanglement parameter ^ for reduced nonlinearity, 
q — Qj/lO. Identical details as in Fig. Q 



2. Entanglement between the two modes 



Writing the quantum state of a TBEC in the form 



JV 



!*(*)> = $> ro (t)|iV>|iv- 

3 

= ^2 c m(*)b'! m ) 



(36) 
(37) 



m=-j 



the well-known von Neumann entropy that shows degree 
of entanglement can be written[l5|: 

E{t) = - * |c m (t)| 2 log 2 |c m (t)| 2 (38) 



where the normalization factor log 2 (iV + 1) _1 was in- 
cluded so that < E < 1 . The expansion coefficients are 
given by 

c m (t) = (m\U(t)\iP(0)) (39) 
= (m\&(t)e- lH,t R(0)\e,^), (40) 

where the rotation operator R, reduced Hamiltonian H' 
and the CSS \6, <j)) are as defined above. Using Eq. 
and inserting the completeness relation X^m' l m ')( m 'l = 
1 one can write 

Cm {t) = ^^mUe-^^-At) 

-i[2m'Slt-qm' 2 t/2] (41) 



where 



xe 



f£L(t) = (m\RHt)\m'} 



(42) 



represents the matrix element of the rotation operator. 
This can be evaluated by first applying the disentangling 
theorem on the rotation operator J?|18|: 



R(X, At) = cxp 



-iAt 



J + e iAt ) 



e Tj+ e wi+\Tr)j, e -T*j- 



(43) 
(44) 



where r = e lAt tan(— -|). Taking into account the unitary 
nature of the rotation operator w (6, (/>) — R{—9,4>) and 
expanding the exponential operators containing J + and 
J- one can write 



T W ,(t) 

• r m,m l y) 



(-r)" r 



L — ✓ 77 I 



Tra + M ) 



1\m! — n' 



indicates a potential for the quantum control of entan- 
glement properties of a TBEC via externally adjustable 
parameters. 



x(to|J"J™ \rri). 



(45) 



The final expression for the matrix element, which was 
obtained by using the ladder operator nature of J± and 
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found to be identical to the irreducible representation of 
full rotation group is: 



A = q. 



F {X) At) 



E 



-iy 



-m-\-n 



U 



n)\ 



(n + m! 



m)\n\ (j 



(j + m)\(j + m'y 



x sm 



(j - m)!(j 

2n^m' —m 



m'y. 

A 
2 



1/2 



i(m— m 1 ) At 



The normalized entanglement parameter i£(f) is plotted 
in Fig. |2| for the four combinations of q and A: q = 0, 
g = qj and A = qj, A = lOgj. It is seen that the non- 
linearity is crucial for achieving high degree of entangle- 
ment. The effect of A was, as expected, to increase the 
frequency of oscillations. To generate Fig. we used 
Eqs. Pfy. igjl , and © with TV = 40 due to the compu- 
tational limitations in numerically calculating factorials 
of large numbers. However, it was found that, with non- 
zero q, as long as the ratio of the various parameters such 
as q/N is kept the same, there are no discernible changes 
in the plot as a function of N . Our result agrees with the 
numerical result obtained by Tonel et aL[28j with higher 
number of atoms. 

It was found that, for the initial phase state \6 = 
7r/2, <j) = 0), the entanglement parameter docs depend 
on the number of atoms, as the normalized coefficients 
|c m (t)| 2 are given in this case by the binomial distribu- 
tion, namely, |c m (t)| 2 = C^ +m . It is also notable that the 
coefficients are no longer time dependent. As the number 
of particles N and hence j increases, the broad binomial 
distribution approximates narrower and narrower Gaus- 
sian as a function of m, and hence the state becomes 
less and less entangled. For example, for N — 400, it is 
found that the scaled E = 0.62. In Fig. we plot the 
entanglement parameter as a function of the number of 
atoms N for the initial phase state. It is seen that the 
value tends towards 0.6 as the number is increased up to 
1000 atoms. Such consistent behavior exhibited for the 
bipartite entanglement suggests a possibility to use this 
property to create a macroscopic matter-wave state with 
known entanglement. 
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FIG. 9: Scaled von Neumann entanglement parameters E for 
N — 40. Top row: q — 0, Bottom row: q — qj. Left column: 
A = qj, Right column: A = 10qj 



FIG. 10: Scaled von Neumann entanglement parameters E 
as a function of the number of atoms iV for the initial phase 
state. 



IV. CONCLUSION 

We have studied the dynamical evolution of a TBEC 
in the presence of an adiabatically varying, off-resonant 
atom-light coupling. The macroscopic, fictitious spin vec- 
tor was found to undergo a rather complex motion, and 
further complications in the form of collapses and revivals 
were found with nonzero nonlinearity. As the main re- 
sult of this paper, we used an exact solution to the time- 
dependent Schrodinger equation to calculate analytically 
the amount of quantum mechanical squeezing and entan- 
glement in a TBEC under various conditions. In partic- 
ular, we considered entanglement generated between the 



atoms and that between the modes. For the case with the 
nonlinearity turned off, although the variance was found 
to be reduced below the SQL (i.e. squeezed), the sys- 
tem never demonstrated atom-atom entanglement, and 
the entanglement between the two modes remained well 
below its maximal value. With the nonlinearity turned 
on, it was found that the atom-atom entanglement can be 
generated initially, and then become unentangled rapidly. 
On the other hand, the entanglement between the two 
modes was found to more or less maintain maximal val- 
ues throughout, except for some fluctuations. The dy- 
namics of entanglement was found to be controllable via 
various parameters present in this system, namely the 
nonlinearity q and the detuning A. 



11 



Potentially useful results of this work are the squeezing 
below the SQL in the atom number difference and relative 
phase which may help reduce the projection noise in spec- 
troscopic and interferometric applications as discussed by 
Wineland et al. 9], and the consistently high degree of 
entanglement between the two modes with the nonlin- 
earity turned on, which could be useful in the context of 
quantum information science. In addition, the dynam- 
ically stable features of the self-trapping state and the 



time-independent entanglement between the two modes 
which is also relatively insensitive to the changes in the 
atom number (for large N) could prove useful in design- 
ing a robust, macroscopic bipartite quantum state with 
known entanglement. Future work could involve study- 
ing quantum control methods for maintaining optimal 
squeezing and entanglement in a TBEC, for such appli- 
cations as precision matter- wave interferometry using en- 
tangled BECsHH. 
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